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Within the framework QCD sum rules, we use the least square fitting method to 
investigate the first radial excitations of the nucleon and light mesons such as p, 
K*, TT , if. The extracted masses of these radial excitations are consistent with the 
experimental data. Especially we find that the decay constant of 7r(1300), which is 
the the first radial excitation of tt, is tiny and strongly suppressed as a consequence 
of chiral symmetry. 
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INTRODUCTION 


The method of QCD sum rules has been widely used to extract the resonance information 
in hadron physics [ij . This formalism is usually applied to study the ground state in a specific 
channel due to the limitation of theoretical accuracy and difficulty of numerical analysis. 
The excitations of mesons have been studied within finite energy sum rules in the literature 
j2-4|. Recently, there are some attempts to study the excitations of the heavy-light mesons 
using the method of QCD sum rules 

The radial excitations have the same spin-parity as the ground state. Experimentally 
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many radial excitations of mesons and baryons have been established (g. Sometimes it is 
quite difficult to identify the radial excitations of hadrons. For example, the situation of 
the radial excitations of the vector charmonium above 4 GeV becomes quite unclear after so 
many charmonium-like XYZ states have been reported experimentally in the past decade. 
Theoretical investigations of the radial excitations are also very challenging. 

In this work, we shall study the first radial excitations of the light mesons and nucleon 
within the framework of the QCD sum rule formalism. We explicitly keep two poles in the 
usual spectrum representation. Then we employ the least square method in the numerical 
analysis to extract the resonance information of the first radial excited state. The extracted 
masses of the radial excitations of the light mesons and nucleon agree with the experimental 
data quite well. 

The paper is organized as follows. In Section [TTl we introduce the QCD sum rule formalism 
and our least-square method. The numerical results are presented in Sections IllltIVIIl The 
last section is a short summary. 


II. FORMALISM 


Within the framework of the QCD sum rule approach, we study the correlation function 
at the quark level 



( 1 ) 


where j {x) is the interpolating current with the same quantum numbers as the hadrons. 
The above correlation function satisfies the dispersion relation 



( 2 ) 


At the quark gluon level, the correlation function can be calculated with the operator product 
expansion. The gluon and quark condensates appear as higher dimensional operators in this 
expansion. At the hadron level, the spectral density of the correlation function can be 
expressed in terms of the hadron masses and couplings. Due to the quark hadron duality, 
we get an equation called the QCD sum rule which relates the correlation function at the 
quark gluon level to the physical states. After making Borel transformation to the sum rule 
in the momentum space, one gets 



( 3 ) 
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where M is the Borel parameter. 

The spectral density usually takes the one-pole approximation 


p('S) — Imll (s) f6 fs Tfl ^ -|- Pcontinuum^ ("S "Sq) ) 
71 


(4) 


where m is the mass of the ground state and sq is the threshold parameter. Above sq, 
the spectral density at the hadron level is replaced by the spectral density derived at the 
quark-gluon level. Now the sum rule reads 


fe-^ = n' (m2) - / e-‘/M^pOPE 

J SQ 


(5) 


The usual numerical method in QCD sum rule analysis is to differentiate Eq. ([5]) with 
respect to and divide the resulting equation by Eq. ([5]) 


° 


J^o g s/m^ gpOPE 

h e-VM^pOPE(5) ds • 


( 6 ) 


One usually plots the variation of the mass versus M^ and sq to find a working window. 

However, the method described above can only be applied to the ground states. In order 
to extract the resonance information of the first radial excitation, we modify the above 
spectral density and explicitly keep the pole of the first radial excitation in the spectrum. 
Now the modified spectral density reads 


p{s) = -ImH (s) = fi6 (s - m2) (s - m'2) pcontinuumd (s - So) • 

TT 


(7) 


To simply the numerical analysis, we use the zero width approximation for both the ground 
state and first radial excitation. The parameters /i and /2 are related to the coupling 
parameters while m and m' are the masses of the ground state and the first radial excitation 
respectively. Now the sum rules read 

/ Pground (s) ds + / e'^/'^Vexcitation (s) ds -h /“ e-'^/'^Vcontinuum (s) ds 

_ j^/ _ j^/perturbation /condensates 

The usual numerical method cannot be applied here because the modified spectrum has 

n 

two mass parameters. We use the least square method |7| to fit these masses and decay 
parameters. The detail of the method are described below. 

As usual in the sum rule analysis, one has to find an optimal working interval of the 
Borel parameter M'^. The lower boundary of M2 is chosen to ensure the convergence of 
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the operator product expansion while the upper boundary is chosen to make the continuum 
contribution remain subleading. 

To get an optimal interval of the Borel parameter we set 


fs7 ® "^'^Vcontinuum (s) ds 


n'(M2 


< tti 


( 8 ) 


which ensures that the continuum contribution remains subleading and determines the upper 
boundary and 

/condensates 

n' (M2) 

which ensures that the OPE is reliable and determines the lower boundary. The two bound¬ 
aries determine the optimal interval of for our numerical analysis. 

The number ai and 0:2 is chosen to ensure a rational contribution of continuum and 
higher order OPE terms. For the meson case, we set ai = 0:2 = « to get a reasonable 
interval of M^. We use different values for ai and 0:2 in the nucleon case. Note that we 
always try a smaller a in the excitation case since the continuum contribution decreases as 
the threshold parameter sq increases. If no reasonable interval of can be got in any way, 
the sum rule may not be appropriate in our numerical method. 

We rewrite the sum rule as 


< 02 


(9) 


_/" 0 ^ Pground ("S) ds -j- J" G ^ Pexcitation ('S) ds 

= g (M^, So) = n' (M^) - pcontinuum («) ds 


which separates the part of expression with physical parameters from the part with just 
Borel parameter and threshold sq. 

With the above expression of g{M‘^,SQ), we can generate a series of points 
{{Mf, g [Mf, So))} by choosing a set {Mf} within the optimal interval of M^. We uni¬ 
formly choose N points in the optimal interval of M^. The number N is chosen to be 20 or 
even larger. 

With the sets {{Mf, g [Mf, sq))}, we use the least-square method which minimizes the 
sum of the squares of the difference between the two sides of the sum rules 

...2 _/2 1 2 


N 

Z 

2=1 


fie ^ + f 2 e 


9 {Mf, So) 


N 


= mm 


( 10 ) 
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to get the best fit of the resonance parameters of the ground state and hrst radial excitation. 

The masses of the ground states of the light mesons and nucleon are measured precisely 
experimentally. The extracted masses from the traditional QCD sum rule formalism with 
the one-pole approximation agree with the experimental data very well. In our analysis 
we hrst use the least square method to reproduce the resonance parameters of the ground 
states. As expected, the resulting masses are consistent with experimental data and those 
extracted from the traditional QCD sum rule analysis. 

Then we use the extracted masses of the ground states as inputs to extract the resonance 
parameters of the radial excited states since the less parameters in the htting will cost less 
computing resource and lead to relatively more stable results. Moreover, we do not £x the 
masses of the ground states in Eq. (1101) in our numerical analysis. Instead we allow them 
to vary around the experimental central value within ±5%. In this way, we extract the 
resonance parameters of the hrst radial excited states numerically. 

We analyze several light mesons and nucleon in the following section. The sum rules 


a 


The nucleon sum rule with the 


of the light mesons can be found in the pioneer paper 

radiative corrections can be found in Ref. jsj. We collect these sum rules in the appendix. 

In our analysis we use the following values for the various condensates and param¬ 
eters [l|, y, 9|: (gg) (2GeV) = — (277tJoMeV)^, (0| |0) = —= 

-1.7 X 10-^GeV^ m, (2GeV) = (95 ± 5)MeV, W>/(g-g) = 0.8 ± 0.3, (0| ^G“^G“J0) = 
0.012+o;o?2GeV‘^, (0| - dyoys^d)^ |0) = fos (0| gg |0)^ ~ 6.5 x 10"^GeV^, 

(0| a* (uy^ysru - dyayst^d) |0) ~ -YW(0|gg|0)^ ~ -6.5 x lO-'^GeV^, 

a, (Q2) = 4Y(Mn(QV)), A = O.lGeV, a, {mz) = 0.1184 ± 0.0007, a, (1.5GeV) = 0.353 ± 
0.006. 


III. THE p MESON 


The interpolating current for the p meson is 




( 11 ) 


and the resulting sum rule can be found in the appendix. The usual single-pole spectral 
density reads 

pG) (s) = 67rVp5 (s - ml) + ^(^ + 0 {s - sq) (12) 
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so[GeV2] 

1.2 

1.3 

1.4 

1.5 

1.6 

MLnlGeV^] 

0.43 

0.43 

0.43 

0.43 

0.43 

M^,,[GeV2] 

0.74 

0.82 

0.88 

0.94 

1.00 

m[GeV] 

0.74 

0.75 

0.75 

0.76 

0.77 

/p[MeV] 

187 

190 

193 

197 

201 

/i[GeV2] 

2.06 

2.13 

2.22 

2.30 

2.39 

min[GeV^] 

10"® 

10“^ 

lo-'^ 

10-5 

10-5 


Table I: The mass and decay constant of the p ground state with a = 0.2 and N = 40. 

We also need the double-pole spectral density 

p^p) (s) = (s - ml) + 67rVp/5 {s - ml,) + ^ (^1 + ^ (s - Sq) , (13) 

where fp and /p/ are dehned as 

(0 \q'y^,q\ p) = mpfpCp, (0 Igy^g] p) = mpfpe^ (14) 

where q = u,d. 

We hrst use the least square method and the traditional one-pole spectrum representation 
with a = 0.2 and = 40 to extract the mass and decay constant of the p meson. The 
results are listed in Table [Tl The parameter fi is related to the decay constant in Eq. (01). 
The values of "min" are the sum of the squares of the differences in Eq. ffTOj) . Only when 
the value of "min" is much smaller than the parameters /f, /| etc, the ht and the extracted 
decay constants are reliable. 

We collect the htting results with the double-pole spectrum in Table [TTl Note that the 
parameter m in Table [TTl is the input to extract the information of the excited state. We use 
a = 0.1 in this case. The threshold sq plays the role of including the hrst radial excitation 
in the spectrum while excluding the contribution from the higher excitations. To check the 
consistency of our htting and dependence of our results on sq, we vary sq in a range. A 
reliable htting requires that the mass m' and decay constant fp of the hrst radial excitation 
should not vary too much with sq. 

From Table [T| we have 


m = (0.76 ± 0.01) GeV, fp = (194 ± 6) MeV 


( 15 ) 
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so[GeV2] 

2.3 

2.4 

2.5 

2.6 

2.7 

M^,JGeV2] 

0.50 

0.50 

0.50 

0.50 

0.50 

M^,,[GeV2] 

1.00 

1.04 

1.10 

1.14 

1.18 

m[GeV] 

0.76 

0.76 

0.76 

0.76 

0.76 

m'[GeV] 

1.24 

1.29 

1.35 

1.38 

1.40 

/p[MeV] 

196 

197 

198 

198 

198 

/p'[MeV] 

130 

141 

152 

161 

170 

/liGeV^l 

2.3 

2.3 

2.3 

2.3 

2.3 

/2[GeV2] 

1.0 

1.2 

1.4 

1.5 

1.7 

min[GeV^] 

10-6 

10-6 

10-6 

10-6 

10-6 


Table II: Masses and decay constants of the p ground state and first radial excitation with a = 0.1 
and N = 40. 

which agrees with the p meson mass from PDG m = 0.77GeV and the experimental 
measurement of the p meson decay constant 0 

yexp _ 216 (5) MeV. (16) 

In order to reduce the dependence on the threshold parameter sq, the extracted values of m 
and fp are the average values of the numerical values in Table [11 From Table jlT] we have 


m' = (1.33 ± 0.07) GeV, fp = (197 ± 1) MeV, fp = (151 ± 16) MeV. (17) 

From PDG, the mass of the first radial excitation is m' = 1.47 GeV and its width is F = 0.40 
GeV. Our extracted p' mass is consistent with the experimental data. At present, the decay 
constant of p' has not been measured yet. 


IV. THE TT AND Ai MESONS 

We adopt the axial current for the pion and Ai mesons 

7^ = (18) 


and the resulting sum rule can be found in the appendix. Besides the Oi pole, the pion also 
contributes to this sum rule due to the partial conservation of the axial vector current. As 









so[GeV2] 

1.3 

1.40 

1.50 

1.60 

1.70 

MLnlGeV^] 

0.52 

0.52 

0.52 

0.52 

0.52 

MLAGeY^] 

1.16 

1.20 

1.28 

1.36 

1.44 

mAi [GeV] 

1.14 

1.18 

1.22 

1.26 

1.28 

MMeV] 

134 

135 

136 

137 

137 

/AjMeV] 

124 

139 

153 

166 

175 

/i[GeV2] 

0057 

0.057 

0.058 

0.059 

0.059 

/2[GeV2] 

0.048 

0.060 

0.074 

0.087 

0.096 

min[GeV^] 

10-^ 

lO"’^ 

10-^ 

10“® 

10“® 


Table III: The mass and decay constant of the Ai meson. We use the least square method and the 
traditional spectrum representation with a = 0.3 and N = 80. 

a Goldstone boson, the pion mass is tiny. Especially in the sum rule analysis, is much 
much less than the Borel parameter M^. We can safely ignore the pion mass and let it be 
zero in the numerical analysis. 

The usual spectrum representation is 

P (s) = (s) + vr/1,5 (s - ^ ^ - so) • (19) 

Our modihed spectrum representation reads 

pW (s) = 7rf^6 (s) + TTf^,6 {s - ml,) + (s - + ^ (^ + ^ ’ (^0) 

where f.„i /ai are dehned as 

(0 \jj,\ vr) = iUp^, (0 |j;| 7i') = iU'p'^, (0 |j;| Ai) = mAj'A^el- (21) 

In the htting, we use the least square method and the traditional spectrum representation 
with a = 0.3 and AI = 80 to extract the Ai mass and decay constant. The results are listed 
in Table IIIII 


In order to extract the resonance parameters of the hrst excitation of the pion meson, 
we employ the modihed spectrum and allow /a^ and mAi to vary around the experimental 
data within ±5%. The numerical results are listed in Table m 
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so[GeV2] 

2.0 

2.1 

2.2 

2.3 

2.4 

2.5 

2.6 

MLnlGeV^] 

0.63 

0.63 

0.63 

0.63 

0.63 

0.63 

0.63 

MLAGeY^] 

1.28 

1.34 

1.38 

1.44 

1.50 

1.56 

1.62 

mAi[GeV] 

1.29 

1.29 

1.29 

1.29 

1.29 

1.29 

1.29 

m(r[GeV] 

1.34 

1.36 

1.31 

1.34 

1.41 

1.43 

1.46 

A [MeV] 

121 

122 

123 

123 

124 

125 

126 

/^jMeV] 

248 

248 

248 

248 

248 

248 

248 

A'[MeV] 

0.2 

0.3 

0.1 

2.3 

0.2 

0.7 

0.1 

/i[GeV2] 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

0.05 

/2[GeV2] 

0.19 

0.19 

0.19 

0.19 

0.19 

0.19 

0.19 

/3[GeV2] 

lO"’^ 

10-^ 

10“® 

10-5 

10-^ 

10-5 

10-® 

min[GeV^] 

10"® 

10-® 

10-5 

10-5 

10-5 

10-5 

10-5 


Table IV: Masses and decay constants of the Ai ground state and the first radial excitation of the 
pion with a = 0.2 and N = 80. 

From Table m we have 

tua, = (1.22 ± 0.06) GeV, U = (135 ± 1) MeV, /a, = (151 ± 20) MeV. (22) 

From PDG we have = 1.23 GeV and F^i^ = 0.40 GeV. We note that the Ai mass from 
the fitting is in rongh agreement with the experimental data. The extracted pion decay 
constant agrees with the experimental data 


/“P = 130MeV. (23) 

However, the extracted Ai decay constant is only half of the experimental data jll| 

= 254 (20) MeV. (24) 

To extract the first radial excitation of the pion meson, we use the experimental data of 
the Ai decay constant as input in the numerical analysis. The results are collected in Table 
lYl We have 


= (1.38 ± 0.06) GeV, U = (123 ± 1) MeV, = (0.6 ± 0.8) MeV. 


(25) 
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The resulting mass of the pion radial excitation agrees with the PDG value very nicely: 

= 1.30 GeV and T^/ = 0.40 GeV Note that the extracted numerical value of is 
not reliable since the parameter /| is even smaller than the "min". In this case, we may get 
a upper bound 

l/sl < \/m^ ~ 0.0032GeV2 (26) 

Accordingly, we get the upper bound for /^/ 

< 0.032GeV (27) 


If the value of /tt' is larger than 0.032 GeV, we should be able to extract its value through 
the least square fitting method. 

In other words, our numerical analysis demonstrates that the decay constant of the pion 
radial excitation tt' is much smaller than the pion decay constant around 130 MeV. This 
interesting fact was also noticed by previous theoretical work including lattice simulations 


mm, 


12 


20| . In fact, the suppression of the tt' decay constant is a consequence of the 


chiral symmetry breaking. In the chiral limit, the decay constants of the pion and its radial 


excitations satisfy the following relation 


2l| 


= 0 , 


(28) 


where (n > 1) is the mass of the pion radial excitation. The pion ground state is 
massless in the chiral limit as a Goldstone boson, hence its decay constant can be large and 
nonzero. For the pion radial excitation, its mass is large and nonzero. Therefore its decay 
constant has to vanish, i.e., 7,^ = 0. 


V. THE K* MESON 

The interpolating current for the K* meson is 

= W7/.S (29) 

and the resulting sum rule can be found in the appendix. The usual single-pole spectral 
density reads 

p (s) = 7r/|-.(5 (s - m^*) + ^ ^ ~ -^o) • (30) 
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so[GeV2] 

1.4 

1.5 

1.6 

1.6 

1.8 

MLnlGeV^] 

0.63 

0.63 

0.63 

0.63 

0.63 

ML.[Gey^] 

1.10 

1.18 

1.28 

1.36 

1.44 

m[GeV] 

0.88 

0.89 

0.90 

0.90 

0.91 

/iC‘[MeV] 

202 

206 

210 

215 

219 

/i[GeV2] 

0.13 

0.13 

0.14 

0.14 

0.15 

min[GeV^] 

10-7 

10-7 

10-7 

10-7 

10-7 


Table V: The mass and decay constant of the K* ground state. We use the least square method 
and the traditional spectrum representation with a = 0.3 and N = 2Q. 


5o[GeV2] 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 

MLnlGeY^] 

0.63 

0.63 

0.63 

0.63 

0.63 

0.63 

ML.[GeY^] 

1.40 

1.46 

1.52 

1.58 

1.64 

1.70 

m[GeV] 

0.89 

0.89 

0.89 

0.89 

0.89 

0.89 

m'[GeV] 

1.22 

1.25 

1.27 

1.33 

1.29 

1.37 

/^‘[MeV] 

200 

201 

202 

200 

207 

207 

/x*'[MeV] 

139 

146 

153 

162 

159 

172 

/i[GeV2] 

0.13 

0.13 

0.13 

0.12 

013 

0.13 

/2[GeV2] 

0.06 

0.07 

0.08 

0.08 

0.08 

0.09 

min[GeV^] 

10-® 

10-® 

10-® 

10-® 

10-® 

10-® 


Table VI: Masses and decay constants of the K* ground state and the first radial excitation with 
a = 0.2 and N = 20. 

Our modified spectrum representation reads 

(s) = nf‘^,5 [s - m\,) + [s - m\,,) + ^ 0 {s - sq) , (31) 

where /k* and fx*' are defined as 

(0 bf'l K‘) = mK./*■.£„ (0 lif l| K-’) = 


(32) 
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The results from the first spectrum representation are listed in Table |V] and those from 
the modified spectrum are listed in Table IVTl From Table IVl we have 

m = (0.89 ± 0.01) GeV, = (210 ± 7) MeV. (33) 

From Table |VT] we have 

m' = (1.28 ± 0.06) GeV, fx* = (203 ± 3) MeV, fx*' = (155 ± 11) MeV, (34) 

where m is an input parameter in Table IVIl The decay constant of the K* was measured to 
be jl0| 

/“f ~ 217MeV. (35) 

From PDG, the mass and width of the K*' are m' = 1.41 GeV, F = 0.232 GeV respectively. 
Glearly our extracted fx* from both fitting agrees with the data. The extracted m' is also 
consistent with data. 


VI. THE ip MESON 

The interpolating current for the ip meson is 

(36) 

and the resulting sum rule can be found in the appendix. The usual spectrum representation 

is 

P (s) = {s + + 6 {s - Sq) . (37) 

We also use the modified spectrum representation 

{s) = ^vr /25 (s - ml) + (s - mj,) + ^ ^ “ ^o) (38) 

where and f^/ are defined as 

(0 |s7^s| ip) = (0 |s7^s| ip') = (39) 

We use the least square method and the traditional spectrum representation with iV = 20. 
Note that there does not exist a working interval of for a = 0.2. So we use a = 0.3 here. 
The results from the first spectrum representation are listed in Table lylTI and those from the 
modified spectrum are listed in Table IVIIIl where m is the input parameter in Table IVIIII 
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so[GeV2] 

1.7 

1.8 

1.9 

2.0 

2.1 

2.2 

M^^JGeV^] 

0.87 

0.87 

0.87 

0.87 

0.87 

0.87 

ML.[GeV2] 

1.66 

1.78 

1.90 

2.00 

2.12 

2.24 

m[GeV] 

1.02 

1.03 

1.03 

1.04 

1.05 

1.06 


217 

221 

226 

231 

236 

240 

/i[GeV2] 

0.016 

0.017 

0.018 

0.019 

0.019 

0.020 

min[GeV^] 

10-9 

10-9 

10-9 

10-9 

10-9 

10-9 


Table VII: The mass and decay constant of the ground state. We use the least square method 
and the traditional spectrum representation with V = 20 and a = 0.3 here. 


so[GeV2] 

3.4 

3.5 

3.6 

3.7 

3.8 

3.9 

4.0 

MLn[Ge\^] 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

1.15 

MLAGeV^] 

2.02 

2.08 

2.16 

2.22 

2.28 

2.36 

2.42 

m[GeV] 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

1.00 

m'[GeV] 

1.45 

1.64 

1.52 

1.55 

1.62 

1.50 

1.51 

UiUeY] 

203 

221 

210 

213 

218 

203 

202 

VlMeV] 

215 

215 

222 

226 

230 

240 

246 

/i[GeV2] 

0.014 

0.017 

0.015 

0.016 

0.017 

0.014 

0.014 

/2[GeV2] 

0.016 

0.016 

0.017 

0.018 

0.018 

0.020 

0.021 

min[GeV^] 

10-10 

10-10 

10-10 

10-10 

10-10 

10-10 

10-9 


Table VIII: Masses and decay constants of the ip ground state and the first radial excitation with 
a = 0.2 and N = 20. 


From PDG, the mass and width of the p> ground state are m = 1.020 GeV and F = 0.004 
GeV while m' = 1.68 GeV, F = 0.20 GeV for hrst radial excitation. The decay constant of 


ground state was measured to be 


10 | 


/“P = 233MeV. (40) 

From Table IVIII we have 


m = (1.04 ± 0.02) GeV, = (229 ± 9) MeV. 


(41) 
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From Table IVIIII we have 

m' = (1.54 ± 0.07) GeV, = (210 ± 8) MeV, /^/ = (228 ± 11) MeV. (42) 

The decay constant of the (p meson from both fittings agrees with the data very well while 
the extracted mass of the first radial excitation is in rough agreement with the data. 


VII. THE NUCLEON 


The interpolating current for the nucleon is 

ri = (43) 

and the resulting sum rule Q| can be found in the appendix. The usual spectrum represen¬ 
tation for the nucleon is 


(s) = (s - + p^^tinuum (s) d {s - Sq) (44) 

where 

Pcontinuum (-S) Imfl (s) 

TT 

= A , _L 

4(27r)^ V 12 TT TT (27r)^ 8 \ vr / 9 tt s 

We also use the modified spectrum representation 



P (-s) = (-S - + (^N'^ (s - m''^) + Pcontinuum (s) 9 {s - Sq) , (46) 

where = 327r"^A^, = 327r"^A^, and Atv is the overlapping amplitude of the interpolating 

current with the nucleon state. 


The results from the first spectrum representation are listed in Table |IX] and those from 
the modified spectrum are listed in Table |Xl To get stable results, we have used the nucleon 
mass and = 2.1 from Table IIXI as input in the numerical analysis of the first radial 
excitation. From PDG, the nucleon mass is m = 0.938 GeV while the mass and width of its 
first radial excitation are m' = 1.44 GeV and F' = 0.300 GeV. From Table HXl we have 


m = (0.93 ± 0.03) GeV. 

From Table |X] we have 

m' = (1.50 ± 0.04) GeV 


(47) 

(48) 


which is in rough agreement with the data. 
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so[GeV2] 

1.80 

1.85 

1.90 

1.95 

2.0 


0.7 

0.7 

0.7 

0.7 

0.7 

ML.[GeV2] 

1.52 

1.54 

1.58 

1.62 

1.64 

m[GeV] 

0.89 

0.91 

0.93 

0.94 

0.96 


1.9 

2.0 

2.1 

2.2 

2.4 

min 

10"^ 

10"^ 

10"^ 

10“^ 

10"^ 


Table IX: The mass of the nucleon ground state with ai = 0.8, ol 2 = 0.4 and iV = 20. 


so[GeV2] 

2.1 

2.15 

2.20 

2.25 

2.30 

2.35 

2.40 

MLnlGeY^] 

0.83 

0.83 

0.83 

0.83 

0.83 

0.83 

0.83 

MLAGeV^] 

1.36 

1.38 

1.40 

1.42 

1.44 

1.48 

1.50 

m[GeV] 

0.929 

0.929 

0.929 

0.929 

0.929 

0.929 

0.929 

m'[GeV] 

1.45 

1.47 

1.48 

1.50 

1.52 

1.53 

1.55 

input 

2.1 

2.1 

2.1 

2.2 

2.2 

2.2 

2.2 

f^N' 

0.67 

0.86 

1.06 

1.28 

1.50 

1.76 

2.00 

min 

10"® 

10“^ 

10-5 

10"^ 

10"® 

10“^ 

10“^ 


Table X: The masses of the nucleon ground state and the hrst radial excitation with ai = 0.7, 
a 2 = 0.3 and N = 20. 


VIII. SUMMARY 

In short summary, we have attempted to extract the masses of the first radial excited 
states of the light mesons and nucleon. In our modified hadronic spectral density, we explic¬ 
itly keep the pole of the first radial excited states together with the ground state. Requiring 
that the operator product expansion converge and the continuum contribution be subleading 
leads to the optimal working interval of the Borel parameter M^. Then a series of "data" 
points (or pseudo-data points) were produced within this working interval of M^. Using the 
usual one-pole spectral density, we can extract the mass of the ground state with the least 
square fitting method, which agrees with the experimental data. Then we use these "data" 
points and the mass of the ground state as input parameters to extract the mass and the 
decay constant of the first radial excited state by the least square method, which are in good 
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agreement with the available data. 

The QCD sum rule method has its inherent accuracy limit due to the various approx¬ 
imations adopted within this framework, such as the truncation of the the OPE series of 
the correlation function, the assumption of the quark-hadron duality, the omission of the 
decay width in the spectral density, the factorization of the four quark condensates and the 
uncertainties of the values of the various condensates etc. In our analysis we only include the 
uncertainty from the fitting using the least square method itself. The least square method 
with the modihed spectrum representation allows us to extract useful information of the 
Erst radial excitations, which depends on the accuracy of the sum rules. It will be very 
interesting to explore whether such a formalism can be applied to the other hadrons. 
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QCD sum rules of the light mesons and nucleon 


For the p meson 


Jdse p{s) = |M^ 


as{M) {o\muUU+madd\o') 


M4 


■ 3 ^- m - 


(^o| Os (M7a75i“«-'i7Q:75i“«^)^|o^ 
M6 


-^TT^ 

9 '' 


Os ( U'yat°‘U+d'1at°‘d Y1 97at“<? 

q=u,d,s 


M6 


(49) 


For the vr meson 


fe “/^V(s)ds = 


1 + 


(M) , l,2(0|^gM^G°.|0) 


+ gTT 


~M^ 


4 

-|7r-^Q(, 


47r® Os (o| 1470:75t“<i<i7c. 75 10) 

M6 

UJat'^U+d'yat'^d 


M6 


(50) 


For the K* meson 


/dse-^/-V^*)(s) = f^ 


1 + 


Os(M) 7r^(0|mu7i'ii+msss|( 


1^2 (o|^G“„G“^|o) r,_3 (o|os(M7c75i“«-s7Q75t“s)^|o) 

' 3 ^ aF 


(51) 


-^7r3 

9 '' 


a a I U-yat‘^U+S'yat°-S Y1 970 *“ I? 

q=u,d,s 


M6 


For the p meson 


/ e ( 5 ) ds = 


-l I aa{M) 6m^{M) , Stt^ (0| ms ss|0) 

7 ' 175 


1^2 _ 4^^3 q,^ (o|gg|o)^ 


+ 37r- 


M4 


M6 


For the nucleon 8 | 


(52) 


40 + 44 + 46 + 48 — /3^e + /3^/e 


(53) 
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where 


io 

A4 (M^,W^) 
Ae 


M^E2 


1 + 


a.. 


a. 


TT 


M*W 


bM^Eo 


4L 


1 - 


Qfc 


TT 


53 , 

— In—- 


13 




^ 3iy2\ _;^2 

1 + ^je T^+Me 


5 1 
6^3 


In 


/i2 


+ s 


'Jp 


'Jp 


a = - (27rf (qq) , b = {2nf {^G^) , Pn = i^n)^ X%, (IGeV) ^ 0.37 
^0 = 1- e-", E2 = l-(^l + x+ e"" 

with X = e {x) = 


In (^ 7 a 2 ) 
In {P/k^) 
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